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ABSTRACT 

Maps of low- inclination nearby galaxies in Sloan Digitized Sky Survey u — g, g — r and r — i colors 
are used to determine whether Ultraluminous X-ray sources (ULXs) are predominantly associated with 
star-forming regions of their host galaxies. An empirical selection criterion is derived from colors of H n 
regions in M 81 and M 101 that differentiates between the young, blue stellar component and the older 
disk and bulge population. This criterion is applied to a sample of 58 galaxies of Hubble type SO and later 
and verified through an application of Fisher's linear discriminant analysis. It is found that 60% (49%) 
of ULXs in optically-bright environments are within regions blueward of their host galaxy's H n regions 
compared to only 27% (0%) of a control sample according to the empirical (Fisher) criterion. This is an 
excess of 3u above the 32% (27%) expected if the ULXs were randomly distributed within their galactic 
hosts. This indicates a ULX preference for young, <;10 Myr, OB associations. However, none of the 
ULX environments have the morphology and optical brightness suggestive of a massive young super star 
cluster though several are in extended or crowded star-forming (blue) environments that may contain 
clusters unresolved by Sloan imaging. Ten of the 12 ULX candidates with estimated X-ray luminosities 
in excess of 3xl0 39 ergs s _1 are equally divided among the group of ULX environments redward of Hn 
regions and the group of optically faint regions. This likely indicates that the brightest ULXs turn on at 
a time somewhat later than typical of Hn regions; say 10—20 Myr after star formation has ended. This 
would be consistent with the onset of an accretion phase as the donor star ascends the giant branch if 
the donor is an <;20 M© star. 

Subject headings: galaxies: general — X-rays: galaxies — X-rays: general 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are defined as 
the most X-ray luminous off-nucleus point-like objects 
in nearby galaxies. Their high luminosities (Lx J> 
10 39 ergs s _1 in the few tenths to ~10 keV range) im- 
ply high accretor masses; hence their interest as a pos- 
sible evolutionary link between stellar-mass black holes 
(Mbh ^1 — 20 M©) known throughout the Local Group 
and super-massive black holes (Mbh ~ 10 6 — 10 9 M©) 
ubiquitous in the nuclei of galaxies. 

Something beyond normal single-star evolution in the 
current epoch is required if ULXs are found to have masses 
much higher than ~100 M w jYungelson et all I200S : 
Ohk ubo et all l2006t IHeger et all . 120031 : IHeger fc Wooslev 



2002). Compact objects with masses >100 M© are often 



referred to as intermediate-mass black holes, or IMBHs, to 
emphasize this non-standard evolutionary path. 

Star formation occurs on a wide range of mass scales, 
even in the strongest starbursts, from the common small 
OB associations up to the rare massive compact super- 
star clusters (SSCs). This hierarchy reflects the mass dis- 
tribution of their originating interstellar clouds and the 
different conditions present at the time of star forma- 
tion. For instance, the most massive clusters form only 
in high-pressure regions such as the high density por- 
tions of dwarf galaxies and galacti c nuclei or the shocked 
regio ns of interacting galaxies ([Elmegreen fc Efremovl . 
1997). Presumably, the most massive stars also form 



only in the most massive clusters (e.g., lElmegreenl 12004 
IWeidner fc Kroupal [2001 . 

Perhaps a similar hierarchical association between ULXs 
and their local environments exists with the lower- 
luminosity, presumably lower-mass, ULXs affiliated with 
small unbound OB associations and only the extremely 
luminous, potentially IMBH, objects coincident with 
SSCs. Certainly, the observed correl ation of ULXs with 
global star for mation rate in sp irals (jSwartz et all 12004 
iGilfanov et~aTl 12004 iLiu et all I2006D implies that many 
bright ULXs are associated with the short-lived popu- 
lation of massive stars. ULXs are also more common 
and more luminous in starburst and interacting/merging 
galaxies than in normal gala xies and they occur at ele- 
vated rates in dwarf galaxies (|Swartz et all l2008f ) . Theo- 
retically, the most promising evolutiona ry path for IMBHs 
(see the review by Ivan der Marell . 12004 for other possibil- 
ities) is merging o f massive stars in the most d ense young 
star clusters (e.g.. iPortegies Zwart et alT . l2004h . 

What is needed is a systematic study of the iocai envi- 
ronments of ULXs to quantify any association ULXs may 
have with particular stellar systems and to quantify the 
properties (e.g., age and mass) of these systems. In this 
way additional constraints on the nature of ULXs can be 
inferred. 

Constraints on the age of some ULXs and of their 
surrounding stellar populations have been obtained from 
high-spatial-resolution optical studies of individual objects 
using Hubble or ground-based 8-m class telescopes. These 
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studies have found several exampl es of late or early B 
stars as eith e r likel y mass dono r s (iLiu et all l2002t [2 004 ; 



Kuntz et all 120051: IPtak et all l2006t iTerashima et al 



2006t iRoberts et all I2008T) or populating; the nearby en- 
vironment (e.g.. ISoria etail 120051 : iGrise et all 12008). In 
some cases, emission nebulae have been found associate d 
with ULXs (|Pakull fc Mirionil . |2002| ; iPakull et all I2006D . 
These results suggest the more common ULX environment 
is a small OB association rather than the compact core of 
a massive SSC. That is, if the hierarchical argument has 
any merit, the observational results to date argue against 
an IMBH association for ULXs. 

However, these studies are feasible only for a few nearby 
ULXs and are limited to those ULXs with luminous donors 
or at least to those located in OB associations. Hence, they 
may not be representative of the average ULX population. 
In fact, there are ULXs in elliptical galaxies. These likely 
have low-mass donors fainter than cur r ent de tection limits 
(e.g.. IPtak etall 120061 : IRoberts et all 120081) . 

Here, a different approach is used to quantify the 
strength of the association between individual ULXs and 
local stellar populations that is applicable to larger sam- 
ples of galaxies and out to larger distances where confusion 
prevents studies of individual stars around a ULX. Specif- 
ically, our goal is to determine what fraction of ULXs 
are associated with Hn regions (currently star- forming 
regions, or with an age <; 10 Myr since the end of star 
formation), and what fraction is spatially coincident with 
young, massive star clusters. To do so, we rely on area- 
averaged optical colors, instead of individual bright stars, 
using the uniform, high-precision, multi-band photometry 
provided by the Sloan Digitized Sky Survey (SDSS). In 
partic ular, the SPSS u band lies blueward of the I? (4000) 
break (|Bruzuall . ll983f ) known to distingu ish young from old 
stella r populations in SDSS data (e.g.. iKauffmann et all 
120031) while u-roru-j col ors can distingu i sh early- from 
late- t ype galaxies in general (|Strateval . 120011 : iBaldrv et all 
l200l . 

A description of the survey sample, details of the re- 
duction of SDSS images to color maps, our definition of 
local environment, and a summary of the X-ray data used 
in this work are presented in §2. An empirical criterion 
to distinguish star-forming regions from non-star-forming 
regions using SDSS colors is developed in §3 along with a 
complem entary analys is using Fisher's linear discriminant 
analysis (Fisher, 1936). Results of applying these criteria 
are given in §4 and the implications for ULXs is discussed 
in §5. 

2. SAMPLE SELECTION 

The galaxies included in this study are a subset of a 
complete sample of 140 galaxies selected for X-ray study 
of their ULX populations. The X-ray sample was cho- 
sen to form a volume-limited set of galaxies that are in 
both the set of all galaxies within 15 Mpc contained in 
the Uppsala Galaxy Catalog with photographic magnitude 
m p < 14.5 mag and the set of all galaxies listed in the In- 
frared Astronomical Satellite (IRAS) catalo gs with a flux 
/fir > 10~ 13 - 3 ergs enr 2 a -1 . The UGC jNilsonl . fl97l 



complete samples: galaxies with angular diameters >1' 
on the first POSS blue prints and galaxies brighter than 
m p — 14.5 in the Zwicky Catalog of Galaxies and Clusters 
of Galaxies. The IRAS fluxes were obtained from 60/im 
and 100/im flux densities list ed in either the IRAS Cata- 
loged Galaxies and Quasars (|Fullmer k, Londsdald. 1989) 
or the IRAS Faint Source Catalog v2.0 (jMoshir et al. . 
fl99l usi ng the formula /ftr/IO" 11 = 3.25560 + 1.265i 00 
following iHo et all (|1997f ) where 5*60 and 5ioo are the total 
flux densities at 60 and 100 fxm, respectively. The IRAS 
catalogs are complete to approximately 1.5 Jy for point- 
like sources. The combined selection criteria favor nearby 
(D < 15 Mpc) predominately optically-bright galaxies 
with at least a modest amount of recent star formation 
as inferred from their FIR luminosities. The sample is 
know n to exclude smaller dwarf galaxies in the neighbor- 
hood (jSwartz et alll2008D . 

Eighty-nine of the 140 galaxies comprising the X- 
ray sample are within the footprint of the 6th 
data release of the S l oan Digitized Sky Survey 
(|Adelman-McCarthv et all . I2008D . Of these, only the 58 
galaxies with inclination i < 65° make up the present 
study as it is too difficult to isolate local features in higher 
inclination galaxies. 

Table 1 lists relevant properties of the galaxies in- 
cluded in this study. Dist ances are ad opted from the 
Nearby Galaxies Catalog (|Tullvl . 119881 ) augmented by 
those computed using the heliocentric velocities listed in 
the NASA/IPAC Extragalactic Database (NED) assum- 



ing H = 73 km s Mpc 



Galaxy celestial coor- 



/FIR dL J-U ergs cm t 
contains all galaxies north of B1950 5 = — 2°30' in two 

4 http: / /www. sdss.org/dr6/tutorials/retrieveFITS. html 
5 http: // montage, ipac . caltech.edu / 



dinates were taken from the NED. Inclinations, optical 
shape (Z?25- major isophotal diameter measured at sur- 
face brightness level of 25 mag s~ 2 in B, position angle, 
and major-to-minor axis ratio) and morphological type 
are from the Third Reference C atalogue of Bright Galax- 
ies (jde Vaucouleurs et all Il99ll ) Table 1 is composed of 
two sections: The 27 low-inclination galaxies hosting one 
or more ULX and the 31 low-inclination galaxies without 
ULXs. 

2.1. The SDSS Data 

All galaxies in the present study are extended objects 
by definition. Therefore, SDSS pipeline products such as 
object lists are not of relevance here. Instead, we begin 
with corrected imaging frames and use these to produce 
color maps for our analysis. Details of this process are 
given in this subsection. 

2.1.1. Preparation of Color Maps 

Corrected frame images of the target galaxies in each 
of the u, g, r and i filters were retrieved from the SDSS 
archival These are the imaging frames with fiat field, 
bias, cosmic ray, and pixel defect corrections applied. Each 
covers an approximately 13.5' x 9.8' region of the sky. In 
several cases, one or more adjacent frames were needed 
to cover a target galaxy; however, each image frame was 
processed separately. 

There are often slight offsets among the 4 filter frames 
of a given field. The Montage(v3.0) image mosaic toolkit^ 
module mProjectPP was used to re-project the frames and 
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module mSubimage was used to center and to crop a small 
border off of each frame. This ensured that the same sky 
coordinates were assigned to a given pixel in all 4 frames of 
a field. Checks against known object positions were made 
to confirm that reprojection did not compromise the as- 
trometric accuracy of the frames. This is important when 
registering ULX positions obtained from X-ray data to the 
final color maps. 

An estimate of the sky level was then made by construct- 
ing an object mask that includes the D25 ellipse of the tar- 
get galaxy, any bright objects elsewhere in the fram e as de - 
termined by our source-finding software (jTennand [2006) , 
and a small border zone to include any artifacts intro- 
duced by the reprojection. The sky level was then defined 
as the average of the unmasked pixels in each frame (after 
properly accounting for the soft bias and converting pixel 
data numbers to true counts using the gain parameters 
provided in the accompanying field information files). 

Sky-subtracted corrected-frame data numbers were then 
converted to calibrated magnitudes using the zero points, 
extinction coefficients and airmass values also provided 
in the field information files. Dereddening to correct for 
Galactic extinc t ion w as also applied using the Hi maps of 
iKalberla etal] (|2005l ) accessed through the FTOOL util- 
ity nh and applying the extinction conversion for each 
SPSS filter as given in the early data release paper 
(|Stoughton et all |2002| ). 

Finally, u — <?, g — r, and r — i color maps were con- 
structed by subtracting the appropriate magnitude image 
in one spectral band from that in another. 

2.1.2. Defining Local Regions 

Local regions are defined here to be the 100x100 pc 2 
area centered on the X-ray position of the ULX. This phys- 
ical size was chosen to be roughly that of a moderate-sized 
H 11 region. For comparison, the aver age diameter f or the 
248 H 11 regions in M 101 identified bv lScowen et all (|l992h 
is 140 pc and the average half-l ight diameter of t he 492 H 11 
regions in M81 compiled by I Pet it et al.1 (|1988f ) is 56 pc. 
Similarly, OB stellar associations in Local Group galax- 
ies are typically 80— 10 pc diameter (e.s., iHunterl . Il995t 
iGouliermis et alll2003D . 

The colors of local regions arc defined as the average 
of the colors of all pixels within the region. Note this is 
not equivalent to the color derived from the difference of 
the average magnitudes in two spectral bands, i.e., it is 
not a flux-weighted but an area-weighted average color. 
In this way, small bright objects are not overemphasized. 
Our averaging method results in colors in units of mag per 
100x100 pc 2 region although, henceforth, color values in 
this work will be quoted simply in magnitudes with the 
fixed area implied. Note this "aperture photometry" also 
differs from that of isolated point-like objects for which a 
local background can be defined using a nearby or sur- 
rounding source-free region. Here, neighboring regions 
contain similar sources of emission and should not be con- 
sidered a background. Thus, only the mean sky values are 
considered background in this work. 

The 100x100 pc 2 region, or resolution element, varies 
from 1.4" x 1.4" (3.6x3.6 pixels subtending 0.396 s pixel -1 ) 
for the most distant galaxy in our sample up to 9.8" x 9.8" 
(24x24 pixels) for the nearest target. In order to maintain 
consistency when comparing local regions hosting ULXs 



to whole galaxies, the same procedure is applied to all 
100x100 pc 2 regions within the D25 area of each galaxy 
(excluding those <3cr above the sky level). We report the 
resulting average galaxy colors which are area-weighted 
colors defined as the sum over the colors of all resolution 
elements within the D25 ellipse of a galaxy, divided by 
the number of elements. Average galaxy colors are listed 
in Table 1. For reference, the number of resolution ele- 
ments range from 334 for the smallest galaxy (UGC 6850) 
to 53262 for the largest (M 101) in our sample. 

2.2. The X-ray Data: ULX Candidates 

The X-ray data for this study includes pointed ob- 
servations obtained using Chandra, XMM-Newton, and 
ROSAT. Th e source-finding and characterization program 
described in iTennantl (|2006l ) was used to locate potential 
(point-like) ULXs. Source positions were compared to the 
NED coordinates of the host galaxy nucleus to cull po- 
tential AGN fr om the ULX list. T hese are the subject of 
another study (|Zhang et all 120091 ). Similarly, comparison 
was made to catalogs of sources to identify and cull known 
background galaxies and foreground stars. Furthermore, 
visual inspection of the SDSS images was made to identify 
potential interlopers not designated elsewhere. (This final 
inspection resulted in the rejection of one foreground star; 
no known background sources were found coincident with 
ULX candidate positions.) Applying the analytical form 
of th e X-ray backgroun d logiV— log5 distribution complied 
bv lMoretti et al.1 (|2003l ) to the D25 areas of each of the 58 
galaxies in our sample, we expect a total of 1.3 background 
sources bright enough to be classified as ULXs. This num- 
ber is small because the typical angular size of our sample 
galaxies is modest and the galaxies are relatively nearby 
so that ULXs sample only the sparsely-populated high-flux 
portion of the background logTV— \ogS distribution. 

Forty-seven ULX candidates were identified in 27 of the 
58 low-inclination galaxies in our sample. Their basic 
properties are listed in Table 2. All ULX candidates were 
observed with Chandra's ACIS detector and these data 
are used to characterize their X-ray porperties. Some of 
the (Chandra) observations are deep enough that the lu- 
minosities of ULX candidates can be accurately estimated 
using spectral fitting (taken from t he literature or per - 
formed here following the methods of lSwartz et all |2004[ ) . 
Other observations are short exposures that accumulated 
only 20—50 events per ULX. For these short Chandra 
exposures, a straight counts-to-flux conversion is applied 
to estimate source luminosities as indicated in Table 2. 
The conversion ass umes a powe r -law s ource spectrum of 
photon index 1.8 (| Swart z et all l2004f ) absorbed by the 
Galac tic hydrogen co l umn density along the line of site 
(from IKalberla et all . [20051 ). The conversion is roughly 
fx ~ lO -11 ^ ergs cm -2 s -1 where fx is the flux and N is 
the observed count rate of events in the 0.5—8.0 keV band. 
The exact conversion depends on the line-of-sight absorp- 
tion. This counts-to-flux conversion was also applied to 
the deep Chandra observations to test the validity of the 
method; the simple estimates agree with published lumi- 
nosities, where available, to within a factor of four which 
is adequate for our present purposes. 

3. STAR-FORMATION COLOR CRITERIA 
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Two methods are employed to differentiate between the 
young stellar component and the older disk and bulge 
population in galaxies based on their SDSS colors. The 
first ( §3.ip is a simple empirical method that selects re- 
gions bluer in all three colors than Hn regions. The sec- 
ond f £|3.2[) is a more rigorous analytical method that finds 
the linear combination of colors which best separate star- 
forming from non-star-forming regions. The two methods 
give very similar results. 

3.1. Empirical Star- Formation Criteria 

Here we establish, empirically, the SDSS colors of re- 
gions known by independent means to be sites of re- 
cent star formation. For nearby galaxies, an appropriate 
and easily-accessible recent star formation indicator is Ha 
emission from Hn regions. This Ha emission comes from 
recombination of gas photoionized by massive O and early 
B stars and therefore indicates ages less than ~10 Myr. 
These stars themselves emit strongly in blue light. We 
therefore expect Hn regions to be bluer than non-star- 
forming regions in the broad SDSS filters. There are, 
however, two factors that could lead to redder colors; dust 
obscuration and the underlying stellar continuum. Thus, 
for instance, a low-luminosity Hn region located within 
the bright bulge of a galaxy can appear more red than a 
similar region located in the disk. 

To develop an emprical star-formation criteria, we ex- 
amine two well-studied large galaxies in our sample; the 
early-type SA(s)ab spiral M81 (NGC 3031) and the late- 
type SAB(rs)cd galaxy M 101 (NGC 5457). The top three 
panels of Figure Q] show the distribution of colors in the 
central bulge, disk, and cataloged Hn regions (whose lo- 
cations and angular sizes are taken from lPetit etaTlll988l ) 
of M 81 in each of the three colors u — g, g — r, and r — i. 
The average of the colors of the Hn regions and of the 
colors of the galaxy as a whole are listed in Table 3 (where 
errors denote one standard deviation). As expected, the 
H II regions are bluer than the bulge and the disk but there 
is considerable scatter which is due in part to obscuration 
by dust. 

We performed a similar analysis for M 101 (bottom pan- 
els of Figure Q] and Table 3) . We find, first of all, that the 
mean colors of the galaxy are bluer than those of M81 
due to the large contribution to the average colors from 
the (red) bulge of M81 in contrast to the (blue) disk- 
dominated late-type galaxy M101. The means of the col- 
or s of H II regions (size s and locations from the tabulation 
of IScowen et all I1992T ) are also bluer than those of M81 
(Table 3). This indicates that the SDSS colors of star- 
forming regions are, indeed, determined in part by the 
underlying stellar light contribution and will vary from 
galaxy to galaxy. 

To summarize, we may formulate our criterion to dis- 
tinguish star-forming regions based on the following em- 
pirical evidence: (1) Star-forming regions are bluer in all 
three SDSS colors, in general, than are non-star-forming 
regions in the same galaxy. (2) Allowance must be made 
for reddening due to dust. (3) Allowance must be made 
for differences among galaxies due to differences in their 
underlying stellar continua. (4) Allowance must also be 

6 http:/ /tmva. sourceforge.net/ 
7 http:/ /root.cern.ch/ 



made for low S/N regions for which colors are poorly con- 
strained. 

Thus, we adopt the following criterion to distinguish 
star-forming from non-star- forming regions. The mean 
colors of the M81 Hn regions define the baseline. These 
values are scaled upward by la to allow for dust redden- 
ing. For each host galaxy, these values are scaled further 
by the difference of the mean colors of the host galaxy to 
that of M81 to allow for variations in the underlying stel- 
lar continua. Symbolically, a region in the i th galaxy is 
flagged as a star-forming region if, for all three colors Cjj, 
j = 1,3, of the region 

C jti < (hj + 8h 3 ) - ((C JMS1 ) - (C jti )) (1) 

where hj is the mean color and Shj is the standard devia- 
tion for Hn regions in M81, (Cj.M8i) is the average color 
of all 100x100 pc 2 regions in M81, and (Cj.i) is the aver- 
age color of regions in the i th galaxy. In addition, regions 
fainter than 3<r above the mean sky level in any of the four 
filters are rejected as low S/N regions under criterion (4) 
above. 

Applying Equation [1] to M 101 results in color criteria 
of C^mioi < 1-07, 0.50, and 0.27 mag for u — g, g — r, 
and r — i, respectively. These are close to the average 
colors (plus Ict) of the known Hn regions of this galaxy 
(Table 3). Table 1 (see § 12.1. 2p includes the mean galaxy 
colors, (Cj,i), and the fraction, /hii, of 100x100 pc 2 re- 
gions classified as bluer than typical H II regions according 
to Equation 1 for each galaxy in the sample. 

3.2. Fisher's Linear Discriminant 

The empirical criterion, Equation 1, treats the 3 SDSS 
colors as independent quantities. The possibility exists 
that some other combination of colors or magnitudes may 
better distinguish star-forming regions from other regions. 
To test this hypothesis, we used the "Toolkit for Multi- 
variate Data Analysis with ROOT" or TMVA0, an add- 
on package to the script interpreter ROOT available from 
CERP>C|. This Toolkit provides several methods of mul- 
tivariate analysis to separate classes of objects. In ap- 
plying the evaluation and test phases of the Toolkit, we 
found that li near di s crimin ant analysis, in the form first 
described by iFisherl (|1936fl , performed very well and we 
adopted it here. The Fisher algorithm finds the linear 
combination of attributes that best separates two or more 
classes of objects. In our case, the attributes are the SDSS 
magnitudes and colors and we wish to distinguish the class 
of star-forming regions from the class of non-star-forming 
regions. The algorithm maximizes the ratio of the dis- 
tance between the classes to the variance within the classes 
thereby guaranteeing maximal separability. The output is 
a set of linear coefficients, dj for the input attributes, Xi 
that gives the Fisher index, F = ao + aix\ + 02X2 + • ■ •, 
such that, in our case, F > would denote a star-forming 
region. 

We began with the known class of Hn regions in M81 
as the training set "source" and all 100x100 pc 2 regions 
in M 81 (that are > 3<r above the sky level) as the train- 
ing set "background". Inspection of the resulting Fisher 
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index indicated that SDSS colors were more important at- 
tributes than magnitudes for separating star-forming from 
non-star- forming regions; which makes sense physically. 
Furthermore, there appeared to be little benefit gained 
by including all three colors with g — r being the least 
important. This naturally suggested using just the two 
attributes, u — g and r — i. Retraining using only these 
two attributes results in a Fisher index 



F = 0.374 - 0.179(u - g) - 0.491(r - i). 



(2) 



Equation 2 did a very good job of separating young star 
forming regions in M81 and provides a criterion compara- 
ble to the empirical method developed in H3.ll 

Applying Equation 2 to M 101 classifies most of the 
galaxy as star-forming because, as before, the mean col- 
ors of the early- type galaxy M81 are redder than those 
of M 101. This is shown graphically in Figure [H The left 
panel shows the distribution of M 101 regions using the H n 
regions of M 81 as the training set (Equation 2). The right 
panel shows the distribution resulting after retraining us- 
ing the Hn regions of M101 itself. Both clearly separate 
the star-forming H II regions from non-star- forming regions 
although the distributions differ in several respects. The 
most important difference is the overall shift or offset from 
the fiducial F — separation criterion to F « 0.1 when us- 
ing the M81 training set. Noting that the average Fisher 
index for M81 (computed using the average colors tabu- 
lated in Table 1 to evaluate Equation 2) is —0.036 and that 
for M 101 is +0.074 suggests that the Fisher index can be 
easily compensated for differences in the underlying galaxy 
colors by applying an offset Aj = ({Fi) — (Fusi)) snch 
that a region in the i th galaxy is flagged as a star-forming 
region if F > Aj. Here, (Fi) is the value of F from Equa- 
tion 2 evaluated using the average u — g and r — i colors 
of the i th galaxy and (Fmsi) is the corresponding value 
for M81; (F M si) = -0.036. This results in A = 0.109 
for M 101. This is wholly equivalent to the scaling by the 
mean galaxy colors that is part of the empirical criterion 
(item 3) leading to Equation 1. 

A comparison of the two methods is given by the frac- 
tion of H II regions actually selected as star- forming. These 
are listed in Table 3 for M81 and M 101. 

Finally, a direct comparison of the two methods is il- 
lustrated in Figure [3] It is a color-color diagram for the 
M81 training data. Each point represents the r — i vs. 
g — r colors of a 100x100 p c 2 region in the ga laxy with the 
known H II regions (from iPetit et"al] . |1988() highlighted. 
The diagonal line traces F = from Equation 2 and the 
horizontal and vertical lines represent two of the three 
empirical cuts from Equation 1 and Table 3. Regions 
to the lower left in the figure are designated blue star- 
forming regions according to these critiria. Note that the 
colors of M81 regions at the red portion of the diagram, 
(u — g,r — i) « (1.8, 0.5), correspond to those of K0 main 
sequence stars, as expected, whereas those in the blue por- 
tion, say (u — g,r — i) « (1.0,0.3), do not correspond to 
colors of a single stellar population; u — g = 1.0 corre- 
sponds to B8 to AO stars whereas r — i = 0.3 is the color 
of GO to G5 stars. This is further evidence of the contribu- 
tion to young star-forming regions of the older underlying 
disk and bulge populations. 

8 available on-line from http://cas.sdss.org/astro/en/tools/chart/chart. 



4. RESULTS 

The two criteria developed in the previous section pro- 
vide similar results. We emphasize results from the em- 
pirical method, W3.ll in the following but also quote values 
from the Fisher linear discriminant when appropriate. 

4.1. Quantitative Considerations 

We can divide the 100x100 pc 2 regions surrounding 
the 47 ULX candidates in our sample into three non- 
overlapping groups: Those rejected as being too faint for 
quantitative analysis of their colors; those meeting the 
star-forming criterion; and those that do not meet the 
star-forming criterion. For brevity, these groups will be 
designated as R for rejected, S for star-forming, and N 
for non-star-forming, respectively; they are listed in Ta- 
ble 2. There are 12 ULX candidates in R regions. Using 
the empirical criterion, there are 21 ULX candidates in 
S regions and 14 in N regions. The results are similar 
for the Fisher criterion except that 6 regions identified as 
S by Equation 1 are classified as non-starforming by the 
Fisher criterion and 2 regions identified as N are classified 
as star-forming. Inspection of the SDSS 3-filter imaged 
did not suggest a reason for this discrepancy. However, 
we note that the two methods agree well in the statistical 
sense. 

As a simple check against the possibility of introducing 
biases into our analysis due to large-scale structures within 
ULX host galaxies, we queried a control sample of regions 
in each galaxy and analyzed them in a manner completely 
analogous to that of regions surrounding the ULX can- 
didates. This sample is defined as the 100x100 pc 2 re- 
gions diametrically opposed to the ULX regions in each 
host galaxy. In this way, at least for well-ordered spirals, 
the control group samples regions similar to those of the 
ULXs on large scales but not on the small scales of individ- 
ual star-forming regions. Using the same group definitions 
as given above for the ULX population, the control sam- 
ple is distributed as 21(21) group R objects, 7(0) group S 
objects, and 19(26) group N objects using the empirical 
(Fisher) criterion. This indicates the ULXs are preferen- 
tially located in star-forming regions and that these regions 
are more luminous than their surroundings. 

We can also estimate the probability of a random sam- 
pling of star-forming regions. We applied the empirical 
star- formation criteria to all 100x100 pc 2 regions in the 
27 galaxies containing ULXs. This ranks 31.8% of all re- 
gions bright enough for analysis as star-forming. Thus if 
ULXs occurred randomly, we would expect only 11 of the 
35 to be classified as star forming regions. The chance of 
finding 21 when only 11 are expected is a 3.1cr deviation. 
On the other hand, finding that 7 of 26 control sample re- 
gions, or 27%, are star-forming is wholly consistent with a 
random distribution of control sample regions. The more 
restrictive Fisher criterion classifies only 23.5% of regions 
as star-forming. The chance of finding 17 ULXs in star- 
forming regions when only 8 are expected is a 3.1<r devia- 
tion while finding of the control sample in star-forming 
regions when 5 are expected is 2.3c low. 

4.1.1. Correlations with Star-forming Regions 
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Fully 60% of those ULX regions with sufficient signal 
(21/35) are in star-forming regions. By definition, these 
regions are as blue or bluer than typical Hn regions of 
their host galaxies and therefore are likely of an age also 
typical of Hn regions which is <;10 Myr; the characteris- 
tic lifetime of the least-massive LyC-producing stars, i.e., 
late O to early B stars of about 15 to 20 M Q . 

4.1.2. Correlations with Host Galaxy Morphology 

Ten of the 27 galaxies hosting ULX candidates in our 
sample are early-type galaxies (SO through Sbc). There is 
no preference for objects in any of the 3 groups, R, S, or 
N, to be associated with a particular galaxy morphological 
type. 

4.1.3. Correlations with ULX Luminosity 

The average X-ray luminosity of group S objects is only 
1.7±0.7 in units of 10 39 ergs s _1 . The average of group 
N is 3.5±3.6 and of group R is 3.9±3.8. These numbers 
reflect the fact that most (74%) of the ULXs in the cur- 
rent sample have observed luminosities estimated to be 
<3xl0 39 ergs s" 1 and that 10 of those 12 more luminous 
than this value are equally divided among groups N and 
R. Thus, the average luminosities in these groups are high 
while the scatter in their luminosities is large. 

4.1.4. Comparison to Mean Galaxy Colors 

A comparison of the color distributions of all 
100x100 pc 2 regions in all ULX host galaxies with the 
color distribution of ULX regions and with the color dis- 
tribution of the control sample is shown in Figure 2] The 
three panels present this comparison for the three colors 
u—g, g — r, and r — i. In each panel, the contributions from 
each galaxy are shifted along the abscissa such that regions 
that exactly meet the color criterion of Equation 1 have a 
color of zero. For clarity of presentation, the histograms 
representing the contributions from all regions (solid his- 
tograms) are normalized to 100. Only resolution elements 
brighter than 3<r above the mean sky level, including those 
containing ULXs and the control sample, are included in 
this census. 

Also shown, in each panel, are best-fitting Gaussian 
functions fitted to each distribution. The best-fit parame- 
ter values are listed in Table 4. 

The ULX u — g color distribution is significantly bluer 
than that of their host galaxies and of the control sample. 
The ULX color distributions in g — r and r — i are simi- 
lar to their host galaxies but still bluer than the control 
sample. The ULX color distributions are broad in all 3 
colors indicating the range of environments in which they 
are found from red galaxy bulges and dust lanes to blue 
young star forming regions. The trends also show, again, 
the strong dependence on the underlying stellar popula- 
tion to the determination of colors of individual regions. 

4.2. Qualitative considerations 

We inspected our color maps and the on-line SDSS im- 
ages and consulted the literature for additional informa- 
tion on the optical properties of the environments of the 
ULXs in our sample. We checked whether some of the 
group N objects may be star-forming regions that are 
heavily reddened by dust either local to the region or su- 
perimposed along the line of sight. Only two of the 14 



group N regions can be confidently associated with strong 
dust lanes based on the morphology of their colors. The X- 
ray luminosities of these two objects are among the highest 
in the sample: 11.1 and 11.3xl0 39 ergs s . Visual inspec- 
tion shows that 5 other group N regions are clearly within 
the bulges of their respective hosts. Six others appear to 
be in disk regions near to but not coincident with blue star 
clusters. These are designated "patchy blue" in Table 2 
along with several of the group R objects. (The remaining 
group N region appears to be a bright star cluster. Its 
colors are inconsistent with those of individual stars and 
of background galaxies. It is designated "star cluster?" in 
Table 2.) Thus, at most only a few of the ULXs in group 
N can possibly be reddened group S objects. We also note 
that galaxies hosting group N regions are not systemati- 
cally stronger (dust-rich) FIR emitters than are group S 
hosts. Neither are they more inclined to our line-of-sight 
than are group S hosts. (Of course, some galaxies host 
both group N and group S objects). 

Visual inspections also showed four ULX regions with 
relatively strong emission in the r filter. Two of these are 
designated group R and two are group N regions with g — r 
colors redder than typical of their host galaxy's Hn re- 
gions. Lacking supporting evidence at this time, we spec- 
ulate that these regions are isolated Hn regions whose 
r-band emission is dominated by Ha with only weak con- 
tributions from underlying starlight. These are denoted 
by "Ha?" in the Comments column of Table 2. 

5. DISCUSSION 

We have devised simple criteria, calibrated on the SDSS 
optical colors of known Hn regions in M81 and M 101, to 
distinguish young, star-forming regions in nearby galaxies. 
We have used these criteria to examine the 100x100 pc 2 
regions around ULX candidates to determine if they are 
associated with the young stellar population. ULX envi- 
ronments are designated as "star-forming" if they are as 
blue as, or bluer than a typical H n region in that galaxy 
(after taking also into account the reddening effects of dust 
and of old underlying stellar populations by scaling by the 
mean colors of that galaxy relative to M 81). A useful fea- 
ture of this method is that colors are distance- independent 
quantities, and are measured on distance-independent re- 
gions. 

For the ULX regions bright enough to permit a quanti- 
tative analysis of their colors, we found that 60% are bluer 
than their host galaxy's H n regions in all three SDSS col- 
ors. This indicates a significant correlation between these 
ULXs and young OB associations. However, none of them 
is associated with a compact super-star-cluster within a 
100x100 pc 2 region. This is at odds with models that 
invoke O-star collision and merge r processes in the col- 
lapsed core of a super-sta r-cluster ()Portegies Zwart et all 
12004 iFreitag et all l2006j) . leading to the formation of an 
IMBH powering the ULX. It is possible that some ULXs 
do indeed form this way but this does not appear to be, 
by any measure, a common formation mechanism. Con- 
versely, there are several ULXs located in the bulges of 
early- type spirals. We interpret them as members of a 
second, low-mass-donor, population in contrast with the 
high-mass systems associated with recent and on-going 
star formation. 

Between these two extremes are several ULXs near 
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patchy blue regions of local star formation. At the same 
time, we found that 10 out of the 12 most luminous ULXs 
in our sample (those with Lx J> 3 x 10 39 ergs s _1 ) are 
located in regions classified as faint or non-starforming. 
Perhaps a couple of these luminous ULX regions are heav- 
ily reddened by dust. Note that the probability of finding 
a ULX outside a star-forming region is 0.5532 (26 of 47 
ULXs are classified as either R or N) therefore the prob- 
ability of getting 10 or more non-S' regions in 12 trials 
selected at random from our sample is 4.4%. 

For a more physical interpretation, we recall that our 
calibrator for star formation is the optical color distribu- 
tion of Hll regions. We speculate that the most lumi- 
nous ULXs (or equivalcntly, the most common phases of 
very high mass transfer) are biased towards early B-type 
donors with an initial mass of i=sl0— 15 M Q and an age 
^10—20 Myr (perhaps at the stage where the B star ex- 
pands to become a blue supergiant). In that case, we ex- 
pect very little residual Ha emission from their surround- 
ings, the O stars having already evolved and died. The 
ULX region in this scenario will appear redder than sur- 
rounding H II regions where star formation is still ongo- 
ing. Individual studies of stellar populations around some 
luminous ULXs are consistent with this explanation: for 
example, an age relO— 20 Myr was inferr ed for the stars 
around NGC4559 X-l dSoria et al |. l2005h. and »2 Myr 
for those around NGC 1313 X-2 (|Grise et all [20081) . 

Another interpretation of the observed separation be- 
tween luminous X-ray sources and nearby st ar-forming re- 
gions (jZezas et all l2002t iKaaret et al.U2004D is that X-ray 
binaries are expelled from their place of birth. This can 
occur through a variety of processes such as an asymmet- 
ric supernova that forms the compact remnant of the bi- 
nary or interactio ns between the bina ry and other stars 
in the cluster (see IKaaret et all 12004 for further discus- 
sion). Typical runaway velocities are of order 10 km s _1 



and typical separations between X-ray binaries and star 
clusters are ~200 p c implying ages of ^20 Myrs. However, 
IKaaret et al.l (|2004f ) also find the cluster-binary separation 
is shorter for more luminous X-ray sources suggesting a 
higher system mass (lower kick velocity) for these objects. 
This is inconsistent with the general pattern described here 
where 10 of the 12 most luminous ULXs are further from 
star-forming regions than are the majority of fainter ULXs. 

The possibility that young, relatively high-mass, sys- 
tems may be located outside well-defined star-forming re- 
gions implies there remains an observational degeneracy 
between the <;20 Myr-old high-mass-donor and the Gyr- 
old low-mass-donor ULX populations. This degeneracy 
is compounded by the reverse scenario; an old low-mass 
system may be located within, or along the linc-of-sight 
towards, a star-forming region. 

However, in a more general sense, we have shown 
that regions of nearby galaxies with differing recent star- 
formation histories can be separated using SDSS colors. 
Hopefully, this method can be applied more generally to 
identify and study correlations between classes of sources 
with stellar populations or substructures in nearby galax- 
ies. 
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Table 1 
Properties of SDSS Galaxies 



Galaxy R.A. Decl. Type D D25 u-g g-r r - i /£ 

(h m s) (dm s) (Mpc) (min) (mag) a (mag) a (mag) a 
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0.34 


0.22 


0.14 


NGC 3738 


11 35 48.8 


+54 31 26 


Im 


4.3 


2.5 


1.13 


0.41 


0.18 


0.27 


NGC 3985 


11 56 42.1 


+48 20 02 


SB(s)m 


8.3 


1.3 


1.00 


0.45 


0.25 


0.47 


NGC 4020 


11 58 56.9 


+30 24 49 


SBd 


8.0 


2.1 


0.93 


0.40 


0.21 


0.31 


NGC 4062 


12 04 03.8 


+31 53 45 


SA(s)c 


9.7 


4.1 


1.25 


0.58 


0.33 


0.33 


NGC 4203 


12 15 05.0 


+33 11 50 


SAB0- 


9.7 


3.4 


1.62 


0.73 


0.39 


0.28 


NGC 4207 


12 15 30.4 


+09 35 06 


S 


8.3 


1.6 


1.36 


0.64 


0.33 


0.37 


NGC 4214 


12 15 39.2 


+36 19 37 


IAB(s)m 


3.5 


8.5 


0.88 


0.34 


0.20 


0.23 


NGC 4245 


12 17 36.8 


+29 36 29 


SB(r)0/a 


9.7 


2.9 


1.53 


0.68 


0.38 


0.36 


NGC 4309 


12 22 12.4 


+07 08 39 


SAB(r)0+ 


11.9 


1.8 


1.37 


0.68 


0.35 


0.31 


NGC 4310 


12 22 26.3 


+29 12 31 


(R')SAB(r)0+ 


9.7 


2.2 


1.33 


0.68 


0.41 


0.34 


NGC 4314 


12 22 32.0 


+29 53 43 


SB(rs)a 


9.7 


4.2 


1.55 


0.71 


0.38 


0.28 


NGC 4370 


12 24 54.9 


+07 26 40 


Sa 


10.7 


1.4 


1.62 


0.74 


0.39 


0.40 


NGC 4395 


12 25 48.9 


+33 32 48 


SA(s)m 


3.6 


13.2 


0.78 


0.27 


0.16 


0.26 


NGC 4414 


12 26 27.1 


+31 13 25 


SA(rs)c 


9.7 


3.6 


1.24 


0.58 


0.33 


0.39 


NGC 4491 


12 30 57.1 


+ 11 29 01 


SB(s)a 


6.8 


1.7 


1.41 


0.65 


0.34 


0.37 


NGC 4509 


12 33 06.8 


+32 05 30 


Sab(pec) 


12.8 


0.9 


0.65 


0.21 


0.05 


0.33 
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Table 2 

Properties of ULX Candidates 



R.A. 


Decl. 


Galaxy 




u-g b 


g-r b 


r-i b 




Group d 


Comments 


(h m s) 


(d m s) 


Host 


(ergs s _1 ) 


(mag) 


(mag) 


(mag) 








02 42 38.89 


-00 00 55.1 


NGC1068 


8.0(11.3) 


1.18 


0.73 


0.31 


11.0 


N 


dust lane 


02 42 39.71 


-00 01 01.4 


NGC1068 


3.0(1.4) 


1.07 


0.41 


0.26 


28.2 


S 




02 42 40.43 


-00 00 52.6 


NGC1068 


1.3(0.6) 


1.60 


0.75 


0.32 


47.2 


N 


bulge 


08 01 48.10 


50 43 54.6 


NGC2500 


6.0(c2f) 


0.88 


0.24 


0.04 


4.0 


N(S) 


patchy blue 


08 01 57.85 


50 43 39.5 


NGC2500 


1.2(c2f) 


0.58 


0.16 


0.26 


14.5 


N 


star cluster? 


08 14 37.02 


49 03 26.6 


NGC2541 


4.4(c2f) 


0.74 


0.24 


0.25 


1.3 


R 


patchy blue 


08 53 33.66 


51 19 29.5 


NGC2681 


1.2(0.9) 


1.34 


0.59 


0.38 


2.2 


R 




08 53 35.75 


51 19 17.3 


NGC2681 


0.7(1.8) 


1.06 


0.65 


0.38 


3.3 


S 




09 22 02.22 


50 58 54.2 


NGC2841 


1.1(0.7) 


1.81 


0.87 


0.47 


13.3 


N 


bulge 


09 55 32.97 


69 00 33.4 


NGC3031 


2.9(3.8) 


1.20 


0.69 


0.35 


13.2 


S(N) 




10 18 12.05 


41 24 20.7 


NGC3184 


1.0(1.7) 


1.88 


0.45 


0.30 


2.2 


R 


patchy blue 


10 18 23.00 


41 27 41.7 


NGC3184 


1.1(2.1) 


0.51 


0.29 


0.06 


3.8 


S 




10 25 06.98 


17 09 47.2 


NGC3239 


1.4(c2f) 


0.66 


0.34 


0.14 


13.9 


S 




10 25 08.20 


17 09 48.3 


NGC3239 


1.9(c2f) 


0.55 


0.36 


0.06 


8.1 


s 




11 20 18.31 


12 59 00.3 


NGC3627 


1.2(c2f) 


0.98 


0.50 


0.27 


13.7 


s 




11 20 20.89 


12 58 46.1 


NGC3627 


4.1(c2f) 


1.77 


0.53 


0.32 


1.4 


R 




11 26 07.33 


43 34 06.3 


NGC3675 


1.7(c2f) 


1.73 


0.71 


0.44 


3.0 


N 


patchy blue 


11 52 37.36 


02 28 07.1 


UGC6850 


1.8(c2f) 


0.41 


-0.06 


-0.34 


13.0 


S 




12 09 22.18 


29 55 59.7 


NGC4136 


2.1(1.7) 


0.35 


0.26 


-0.07 


4.5 


S 




12 10 33.76 


30 23 58.0 


NGC4150 


3.9(c2f) 


1.61 


0.75 


0.39 


17.7 


N 


bulge 


12 15 10.91 


20 39 12.4 


NGC4204 


1.6(c2f) 


-0.06 


0.33 


-0.44 


1.8 


R 


isolated H n? 


12 28 17.83 


44 06 33.9 


NGC4449 


0.6(0.8) 


0.17 


0.25 


-0.19 


63.7 


S 




12 30 29.55 


41 39 27.6 


NGC4490 


1.1(1.0) 


0.81 


0.75 


-0.30 


15.5 


N(S) 


patchy blue, Hn? 


12 30 30.82 


41 39 11.5 


NGC4490 


2.7(2.6) 


0.54 


0.30 


0.18 


18.7 


N 


patchy blue 


12 30 32.27 


41 39 18.1 


NGC4490 


2.0(1.7) 


0.74 


0.40 


0.10 


14.1 


S 




12 30 36.32 


41 38 37.8 


NGC4490 


2.1(1.8) 


1.12 


0.37 


0.22 


45.3 


N 


bulge 


12 30 43.26 


41 38 18.4 


NGC4490 


3.8(3.1) 


0.68 


0.26 


0.10 


14.8 


S 




12 35 51.71 


27 56 04.1 


NGC4559 


22.9(9.4) 


0.28 


-0.08 


-0.06 


1.8 


R 


patchy blue 


12 35 57.79 


27 58 07.4 


NGC4559 


2.8(1.8) 


0.76 


0.32 


0.22 


21.4 


S(N) 




12 35 58.56 


27 57 41.9 


NGC4559 


15.0(11.1) 


1.11 


0.59 


0.31 


15.5 


N 


dust lane 


12 36 08.90 


19 19 55.9 


NGC4561 


3.0(c2f) 


0.57 


0.35 


-0.04 


1.9 


R 


patchy blue 


12 41 29.14 


41 07 57.7 


NGC4618 


1.7(c2f) 


0.54 


0.23 


0.09 


4.6 


S 




12 41 52.72 


41 16 31.7 


NGC4625 


1.3(c2f) 


1.00 


0.51 


0.21 


20.5 


S(N) 




12 50 25.70 


25 31 29.8 


NGC4725 


1.6(1.8) 


0.50 


0.60 


-0.10 


9.1 


S 


Hn? 


12 50 26.37 


25 33 19.4 


NGC4725 


2.1(2.0) 


-0.21 


0.41 


1.11 


0.7 


R 




12 50 27.39 


25 30 26.5 


NGC4725 


1.0(0.8) 


2.04 


0.88 


0.48 


6.4 


N 


bulge 


12 50 36.88 


25 30 28.4 


NGC4725 


1.8(1.6) 


1.23 


0.79 


0.52 


0.9 


R 




19 ^ 19 31 


UU U / OL.ij 




9 A(oOf\ 


n OA 


u.uo 


u.uo 


u.u 


o 




13 15 19.54 


42 03 02.3 


NGC5055 


6.0(13.5) 


0.73 


0.40 


0.30 


2.2 


R 




13 15 39.33 


42 01 53.4 


NGC5055 


0.4(0.6) 


1.17 


0.69 


0.34 


8.2 


S(N) 




13 16 02.27 


42 01 53.6 


NGC5055 


1.5(1.3) 


1.12 


0.69 


0.18 


6.5 


S 




13 29 38.61 


58 25 05.6 


NGC5204 


4.7(1.2) 


0.69 


0.23 


0.14 


26.8 


S(N) 




14 03 32.39 


54 21 02.9 


NGC5457 


3.4(2.0) 


0.96 


0.34 


0.23 


8.6 


S(N) 




14 04 14.29 


54 26 03.8 


NGC5457 


1.7(1.2) 


0.52 


0.08 


0.74 


1.2 


R 




14 04 59.74 


53 38 08.9 


NGC5474 


5.1 (c2f) 


0.92 


0.30 


0.25 


8.4 


N 


patchy blue 


14 19 39.39 


56 41 37.8 


NGC5585 


2.8(c2f) 


-0.31 


0.78 


-0.60 


1.1 


R 


isolated H n? 


14 19 52.14 


56 44 17.3 


NGC5585 


1.0(c2f) 


1.08 


0.38 


0.28 


3.7 


N 


patchy blue 



a Observed luminosity estimated from count rates, values from spectral fits in parentheses. "c2f" denotes too few counts 
for spectral fitting. 

b Average colors within 100x100 pc 2 region. 
C SDSS detection significance above sky level in a. 

d Group designation R for rejected, S for star-forming, and TV for non-star-forming according to the empirical criteria, 
^3.11 values in parentheses according to the Fisher criterion, H3.21 when different. 
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TABLE 3 
SDSS Color Distributions 



M81 


u — g 


g — r 


r — i 


f Empirical 


f Fisher 


Galaxy ((Cj, M8 i)) 


1.30±0.34 


0.63±0.11 


0.36±0.09 






Hn (hj±8h 3 ) 


1.03±0.30 


0.61±0.10 


0.31±0.09 


77.0 


64.4 


Disk 


1.26±0.26 


0.62±0.08 


0.36±0.09 






Bulge 


1.66±0.10 


0.74±0.04 


0.41±0.03 






M101 


u — g 


g — r 


r — i 


f Empirical 


f Fisher 


Galaxy 


1.04±0.24 


0.42±0.10 


0.23±0.06 






Hn 


0.93±0.26 


0.42±0.13 


0.17±0.10 


68.6 


69.5 



/Empirical and f Fisher are the fractions of regions classified as bluer than Hn regions 
according to the criteria of i )3.1l and £ 13.21 respectively 



TABLE 4 

Gaussian Fits to Color Distributions 







Galaxies 


ULXs 


Control 






Center 


a 


Center 


CT 


Center a 


u 


-9 


0.040 


0.268 


-0.084 


0.284 


0.135 0.172 


9 


— r 


-0.058 


0.091 


-0.047 


0.145 


-0.013 0.087 


r 


— i 


-0.011 


0.060 


-0.003 


0.069 


0.020 0.045 
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Fig. 1. — Color distributions in the bulge, disk, and H II regions of M 81 (top) and the disk and H II regions of M 101 (bottom). Left-to-right: 
u ~~ 9> 9 — r i an d f — i colors of all 100x100 pc 2 regions in the color map; Heavy line: Hll regions, light line: disk, dotted line: bulge. 




0.2 

Fisher Index 




0.2 

Fisher Index 



Fig. 2. — Left: Image of the central field of M81 with all regions classified as star-forming according to the Fisher linear discriminant 
(F > 0, Equation 2) highlighted. This result can be be compared to the empirical method displayed in the rightmost panel of Figure [T] 
Middle: Distribution of all 100x100 pc 2 regions of M 101 in the Fisher index (heavy line) using the M81 Hll regions as the training set 
(Equation 2). The light line traces the histogram of the known Hll regions in M101. The area under each histogram is normalized to unity. 
Right: Same as middle panel but with the Fisher index trained using the H II regions of M 101. Note the fiducial F = separator applies here 
but that F m 0.1 would be more appropriate for the result shown in the middle panel. See text for further discussion. 
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Fig. 3. — Color-color diagram of the 100x100 pc 2 regions in M81. Known Hn regions (Petit ct al., 1988) are indicated by crosses. The 
horizontal and vertical lines are the empirical delineations (Equation 1) and the diagonal line is the Fisher index, F = (Equation 2) 
separating blue star-forming regions below and to the left of these lines and red bulge and disk regions above and to the right. 




Fig. 4. — Left-to-right: u — g, g — r, and ? — i color distributions for all resolution elements in all galaxies (light lines), all ULX regions (heavy 
lines), and all control regions (dotted lines) with significance >3<j above the local sky background. Smooth curves are Gaussian function fits 
to the histograms. 



